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Two design methods for multi -parameter dynamic systems 
are proposed. They are intended to eliminate the limita- 
tions and disadvantages of the existing design methods. The 
powerful mathematical tools of optimal control theory are 
applied to the practical design problems of classical con- 
trol . 

The first method is intended for linear systems only; 
the design problem is solved in the s-domain, by finding 
"the best root locations" of the system's characteristic e- 
quation. In the second method, the design problem is solved 
by finding "the best response" of the system in the time 
domain. The second method is applicable to a wide range of 
dynamic systems; it can be used to synthesize linear, non- 
linear and sampled-data systems, and systems with time de- 
lay. This method is also extended to a numerical stability 
analysis procedure. 

Fourteen examples are presented to illustrate the 
applications of the methods. 



2 



t POSTGRADUATE SGHUOi- 

EREV, CALIF. 93940 



TABLE OF CONTENTS 



Page 

I . INTRODUCTION 7 

II. LINEAR SYSTEM DESIGN BY OPTIMIZATION FOR 10 

THE BEST ROOT LOCATIONS IN THE s -DOMAIN 



A. 


GENERAL ‘ 


10 


B. 


PHILOSOPHY OF OPTIMIZATION 


13 


c. 


CONSTRAINTS 


15 


D. 


FACTORING THE CHARACTERISTIC EQUATION 


15 


E. 


THE PERFORMANCE INDEX 


18 


F. 


PROGRAMMING 


20 


G . 


EXAMPLE 1 


21 


H. 


EXAMPLE 2 


33 


I. 


EXAMPLE 3 


38 


J. 


EXAMPLE 4 


40 


K. 


CONCLUDING REMARKS 


45 


OPTIMIZATION FOR THE BEST RESPONSE IN THE 


48 



TIME DOMAIN — APPLICATIONS TO LINEAR 



SYSTEM DESIGN PROBLEMS 

A. GENERAL 48 

B. SELECTING THE MODEL RESPONSE 49 

C . CONSTRAINTS ' 50 

D. DEVELOPMENT OF THE METHOD 51 

E. EXAMPLE 5 54 

F. EXAMPLE 6 58 

G . EXAMPLE 7 59 



3 



Page 

H. EXAMPLE 8 65 

I. CONCLUDING REMARKS 68 

IV. OPTIMIZATION FOR THE BEST RESPONSE IN THE 72 

TIME DOMAIN — APPLICATIONS TO NONLINEAR 

AND SAMPLED-DATA SYSTEMS AND SYSTEMS WITH 
TIME DELAY 

A. GENERAL 72 

B. EXAMPLE 9 72 

C. EXAMPLE 10 84 

D. EXAMPLE 11 89 

E. CONCLUDING REMARKS 104 

V. A TIME-DOMAIN STABILITY ANALYSIS METHOD 105 

USING NUMERICAL TECHNIQUES 

A. GENERAL 105 

B. GENERAL PHILOSOPHY 105 

C. EXAMPLE 12 108 

D. EXAMPLE 13 108 

E. EXAMPLE 14 110 

F. CONCLUDING REMARKS 119 

VI. CONCLUSIONS 122 

COMPUTER PROGRAM I 123 

COMPUTER PROGRAM II 128 

BIBLIOGRAPHY 135 

INITIAL DISTRIBUTION LIST 136 

FORM DD 1473 139 



4 



ACKNOWLEDGEMENT 



The author acknowledges with gratitude the advice and guidance 
of Professors George J. Thaler and Donald E. Kirk. 



5 



I. 



INTRODUCTION 



Classical control theory and modern control theory 
have been extended in diverging routes since the inception 
of modern control theory. Modern techniques have been suc- 
cessful for larger scale problems and for problems which 
include some new performance criteria, such as expenditure 
of control energy, etc., but a general purpose or a special 
purpose computer in the control loop is usually necessary 
to implement the solution. 

Most design techniques of classical control theory 
were developed before the extensive application of digital 
computers to system design problems. Using a digital com- 
puter in the design phase of a control system is usually 
feasible, especially if the design can be implemented by 
using passive circuit elements rather than complicated and 
expensive memory units. 

During the last decade mathematical methods and numer- 
ical teqhniques have been developed and applied to the prob- 
lems of "Optimal Control Theory". The same methods have 
found application in other areas of science, but they have 
not been applied efficiently to the design procedures of 
the dynamic systems which have been defined in the sense of 
classical control theory. The classical techniques, al- 
though having many drawbacks and limitations (see Section 
II. A), were well developed by the time effective optimiza- 
tion methods arrived on the scene. The classical methods 
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have been successful for a restricted class of design prob- 
lems; for this reason probably, modern control theory grew 
independently, rather than improving and replacing the 
existing methods. 

In this thesis it is intended to narrow the gap be- 
tween the classical and modern control theories. Powerful 
mathematical tools of optimization theory are applied to 
the practical design problems of classical control. Two 
general design methods for multi-parameter dynamic systems 
are proposed: 

1. Optimization for the Best Root Location in the 
s -Domain, and 

2. Optimization for the Best Response in the Time 
Domain. 

Chapter II of the thesis is devoted to the first meth- 
od which is intended for linear systems only — the design 
problem is solved in the s-domain by finding "the best root 
locations" of the system's characteristic equation. Four 
linear system design examples are included. The rest of 
the thesis is devoted to the second method and its various 
applications. In Chapter III, a general description and 
the philosophy of the method "Optimization for the Best Re- 
sponse in the Time-Domain" is given, and it is applied to 
the same linear system design problems discussed in Chapter 
II, to provide a basis for comparison of the results ob- 
tained by both methods. 
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The time-domain method is applicable to a large varie- 
ty of system design problems. To demonstrate this point. 
Chapter IV is devoted to various applications of the method. 
Three examples for nonlinear and sampled-data systems and 
a system with transport lag are presented in Chapter IV. 

In Chapter V, a numerical method to investigate the 
stability limit of the dynamic systems, is introduced as an 
extension of the time-domain method. 

In the year 1964, the gap between control theory and 
the practical control problems became substantial and ex- 
tensive research effort has been directed to bridge this 
gap between the theory and paractice (see, for example. Ref. 
2) . The methods presented in this thesis can be thought of 
as another approach to suboptimal system design problems. 

In this respect, these methods can be considered as applica- 
tions and extensions of Optimal Control Theory; on the other 
hand, these methods provide better solutions to the practi- 
cal dynamic system synthesis problems of Classical Control 
and they can be considered as the extensions of Classical 
Control. Since the extensions from Optimal Control Theory 
and Classical Control are toward each other, these methods 
hopefully provide one of the necessary bridges to reduce 
the gap . 
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II. LINEAR SYSTEM DESIGN BY OPTIMIZATION FOR 



THE BEST ROOT LOCATIONS IN THE S- DOMAIN 
A . GENERAL 

Present design methods of classical control theory are 
limited and have some serious disadvantages. Root-locus 
and frequency-domain techniques, for example, can handle 
systems which have only one parameter. Algebraic methods 
(in the coefficient and parameter planes) are two-variable 
procedures. The major design tools of these methods are a 
set(s) of graphs, and obtaining these graphs for higher- 
order systems is quite a laborious task. The trial-and- 
error nature of the classical design procedures is, however, 
the most significant limitation of these methods. 

Since the dynamic behavior of second-order systems has 
been thoroughly studied and is well known, most of the pre- 
sent methods for high-order system design depend on the 
dominant-roots concepts. Finding a pair of complex dominant 
roots for a high-order system is equivalent to approximat- 
ing the system by a second-order system. 

A pair of complex roots of high-order system can be 
made dominant by satisfying at least one of the two condi- 
tions given below: 

1. Transients in the time response due to undesired 
roots (all the roots except the selected "dominant" pair) 
last for a much shorter time than the transients of the se- 
lected pair. This is achieved if the time constants of the 
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undesired roots are much smaller than those of the selected 
pair. In the s-domain, the magnitudes of the real parts of 

the undesired roots must be much larger than the magnitu- 

; 

des of the real parts of the selected pair to satisfy this 
condition. The remote undesired roots may cause undesir- 
able effects in the time response for a short period of 
time (near zero time) , but thereafter the time response is 
governed by the selected roots. 

2. The coefficients in the time response associated 
with the desired roots are much larger than those associ- 
ated with the undesired roots. In other words, the contri- 
bution of the selected pair to the transient response is 
much larger than the contribution of the other roots. In 
the s-domain, this condition is satisfied if the residues 
of the desired roots are much larger than the residues of 
the other roots. 

If both of these conditions are satisfied, the time 
response of the selected pair dominates the overall time 
response, and the selected pair of roots is called 
"dominant". If only one of the conditions is satisfied an 
acceptable transient response may still be obtained. The 
first condition alone, for example, has been used in the 
s-domain design procedures to avoid the lengthy and time- 
consuming residue calculations at each step of the trial- 
and-error procedures. 

When designing a high-order feedback system by rbot- 
locus methods, or by algebraic methods, the designer 
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usually selects a pair of complex root locations in the s- 
domain, and then tries to locate two of the roots of the 
system's characteristic equation at these selected loca- 
tions and the other (undesired) roots at remote points of 
the left-hand side of the s-plane. If this can be achieved, 
the transient response of the system closely follows the 
selected (desired) second-order system's response. There 
are usually other specifications to be met, such as the 
steady-state error, maximum allowable gain, and limitations 
on the system components, which are treated as additional 
constraints . 

The number of free independent parameters (N c ) cannot 
exceed the order of the characteristic equation (n) , and if 
the number of free independent parameters is equal to the 
order of the system (N c = n) , the roots of the characteris- 
tic equation can be located at any point in the s-plane. 

The design is completed by simple algebraic manipulations 
(see, for example. Ref. 12). If the number of free indepen- 
dent parameters is less than the order of the characteris- 
tic equation of the system (N c < n) , however, the designer 
loses his control over some of the roots. Since one free 
parameter exactly locates one root, when all parameters are 
used to locate an equal number of roots, the excess (n-N c ) 
root(s) may be anywhere in the s-plane. Some of these 
roots may move very close to the origin and become dominant, 
or some may even go into the right half of the s-plane and 
make the system unstable. 



12 



Two parameters can be used to locate a pair of complex 
roots exactly at the specified locations, but these roots 
can be constrained to move on a line by using just one pa- 
rameter. For example, their real parts (a) or imaginary 
parts (to) can be fixed, or they can be forced to move on a 
constant damping ratio (C)-line or on a constant natural 
frequency ( tO n ) -circle . 

When designing with a fewer number of parameters than 
the degree of the characteristic equation, the procedure is 
a trial-and-error method and success may depend primarily 
on the past experience of the designer. There are several 
methods for locating two of the roots of the characteristic 
equation at the desired locations, but if the number of 
free parameters is less than the order of the system, none 
of these methods ensure 

a, the stability of the system, and 

b. the dominance of the selected roots. 

B. PHILOSOPHY OF OPTIMIZATION 

Since placing a pair of complex roots at the desired 
locations has no significance without ensuring their domi- 
nance and the system's stability, the method proposed in 
the first part of the thesis provides these two conditions 
by confining the undesired roots to appropriate region (s) 
of the s-plane, and then finds the optimum parameter values 
to locate the dominant roots as close to the specified 
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locations as possible^. Since the undesired roots are not 
placed at specific locations (they are only confined to 
stay in a remote region of the s-plane) none of the free 
parameters are "used up"; only their freedom is limited. 

The problem is solved in the "Parameter Space" which 
is a N^-dimensional Euclidean space. A performance index 
(J) , which is usually a measure of closeness of the dominant 
roots to the desired root locations, is defined. By mini- 
mizing the performance index with respect to the selected 
free parameters (z) (see Section II. D) one obtains the 
"best possible solution" to the design problem as it is 
specified. The design specifications yield some linear and 
nonlinear constraining equations and/or inequalities; it 
is assumed that these constraining equations define a bound- 
ed region in the N c ~dimensional parameter space. Some of 
the constraining equations are generally nonlinear and it 
may sometimes be impossible to find a feasible point ( a 
point which simultaneously satisfies all of the constrain- 
ing equations and inequalities) to minimize the performance 
index. In practical terms, this means that the design (or 
compensation) is not always attainable by using a certain 



For the simplicity of explanation it is assumed here 
that the design is carried out for a selected pair of 
complex roots; but this is not mandatory. The method is 
generalized (in Section II. B) and an example is given 
(Example 3) for the design of systems without a pair of 
selected dominant roots. 



type of compensator. If this occurs, the designer should 
change the system configuration (usually by changing the 
type of compensator, or by inserting additional compensa- 
tion) . 

C. CONSTRAINTS 

Performance specifications and available system com- 
ponents dictate the constraints, which can be divided into 
three general groups : 

1. Constraints due to physical, limitations, such as 
limitations on the available system components. 

2. Constraints due to static (or steady-state) per- 
formance specifications. 

3. Constraints due to dynamic performance specifica- 
tions . 

The constraints in the first two groups are linear or 
can be approximated by linear inequalities. The constraints 
in the last group are obtained by factoring the system's 
characteristic equation (as explained in the next section) 
and they are expressed as "n" linear and/or nonlinear si- 
multaneous algebraic equations. 

D. FACTORING THE CHARACTERISTIC EQUATION 

The denominator of the system's transmission function 
(from input to output) is called "the characteristic poly- 
nomial", and when it is equated to zero it is called "the 
characteristic equation". The characteristic polynomial of 
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